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We present results from a hybrid description of Cu+Cu collisions using 3+1 dimensional hydrody- 
namics (music) for the bulk evolution and a Monte-Carlo simulation (martini) for the evolution of 
high momentum partons in the hydrodynamical background. We explore the limits of this descrip- 
tion by going to small system sizes and determine the dependence on different fractions of wounded 
nucleon and binary collisions scaling of the initial energy density. We find that Cu+Cu collisions 
are well described by the hybrid description at least up to 20% central collisions. 
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I. INTRODUCTION 

In order to learn about the range of applicability 
of both hydrodynamics and energy loss mechanisms in 
heavy-ion collisions, wc study Cu+Cu collisions at y/s = 
200 GeV and different centralities. Cu+Cu collisions are 
ideally suited for this because they are systems with a 
number of participants that corresponds to that in pe- 
ripheral Au+ Au collisions [l[ , and can be studied experi- 
mentally with reduced uncertainties Q . One should thus 
be able to study the transition from "soft" (hydro-like) 
features to "hard" features (jet- like). 

Hydrodynamics assumes that the system is equili- 
brated within a short time of the order of tq = 0.5fm/c. 
This assumption could be fulfilled for central Cu+Cu col- 
lisions but is unlikely to hold for peripheral ones. Hence, 
it is important to study a wide range of centrality classes 
to see when and for which px range the assumption 
breaks down. 

In this work we are able to describe the entire exper- 
imentally covered pt range. The soft or low-p^ sector 
of the matter produced is described using 3+1D rela- 
tivistic hydrodynamics (music [H) and the hard sector 
Opt ^ 2 GeV) with the Monte-Carlo simulation MARTINI 
[J], which uses the hydrodynamic calculation as input 
and evolves the hard partons in the soft hydrodynamic 
background with the McGill- AMY (Arnold, Moore and 
Yaffe) energy loss approach which treats radiative 0, Q 
and elastic processes Q- 

We can therefore achieve a complete description of 
Cu+Cu collisions under the assumption that both hydro- 
dynamics and the energy loss formalism are valid. Know- 
ing that we reach the limits of applicability of both mod- 
els particularly for large centralities, this study allows to 
pin down observables that are sensitive to a failure of our 
approach and draw conclusions about the nature of the 
produced system in Cu+Cu collisions. 

II. SOFT PART - HYDRODYNAMICS 

To describe the soft part of the matter produced in a 
heavy-ion collision we employ hydrodynamics. More pre- 
cisely, we use the ideal 3+1D relativistic implementation 
of the Kurganov-Tadmor scheme [|[ - MUSIC - introduced 
in @. 



The hydrodynamic equations have the following gen- 
eral form: 



dtp a 



-V-J„ 



(1) 



where a runs from to 4, labelling the energy, 3 compo- 
nents of the momentum and the net baryon density. 

The task is to solve these equations together with the 
equation of state. The Kurganov-Tadmor (KT) method, 
combined with a suitable flux limitcr, is a non-oscillatory 
and simple central difference scheme with a small artifi- 
cial viscosity that can also handle shocks very well It is 
Riemann-solver free and hence does not require calculat- 
ing the local characteristics. 

When describing heavy-ion collisions it is useful to em- 
ploy t — Tj s coordinates defined by 



t = t cosh rj s , 
z = t sinh r) s . 

Then the conservation equation cLJ^ — becomes 



where 



d T (rJ T ) + d Vs J Vs + O v {tJ v ) = . 



J T = (cosh ri s J° — sinh?7 s J 3 ) , 
J Vs = (cosh rj s J 3 - sinh 77 s J°) , 



(2) 



(3) 



(4) 
(5) 



which is nothing but a Lorentz boost with the space-time 
rapidity 77., = tanh~ 1 (z//;). The indices v and w refer to 
the transverse x,y coordinates which are not affected by 
the boost. Applying the same transformation to both 
indices of T Ml/ , one obtains 



d T (rT TT ) + d rh (T"" T ) + 8 v (tT vt ) + = , (6) 



and 



d T (rT"' s ) + d rjs (T r,B71s ) + d v (tT v ^ ) + T Tr,s = , (7) 



and 



dr{rT TV ) + d Va (T r >° v ) + d w (rT wv ) = . 



(8) 



These 5 equations, namely Eq. ([3]) for the net baryon cur- 
rent, and Eqs. (J6j [3 [8j) for the energy and momentum are 
the equations we solve with the KT scheme described in 
detail in 
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To close the set of equations (0 [6J [3 [8]) we must 
provide a nuclear equation of state V(e,p) which re- 
lates the local thermodynamic quantities. The presented 
calculations all employ an equation of state extracted 
from recent lattice QCD calculations Q. There, several 
parametrizations of the equation of state which interpo- 
late between the lattice data at high temperature and 
a hadron-resonance gas in the low temperature region 
were constructed. We adopt the parametrization "s95p- 
vl" (and call it EOS-L in the following), where the fit to 
the lattice data was done above T = 250 MeV, and the 
entropy density was constrained at T = 800 MeV to be 
95% of the Stefan-Boltzmann value. Furthermore, one 
"data-point" was added to the fit to make the peak in 
the trace anomaly higher. See Ref. Q for more details 
on this parametrization of the nuclear equation of state. 

The initialization of the e nerg y density is done using 
the Glauber model (see Ref. [10[ and references therein) : 
Before the collision the density distribution of the two 
nuclei is described by a Woods-Saxon parametrization 



Pa{t) 



Pa 



l + exp[(r-i?)/d] 



(9) 



with R = 4.163 fm and d = 0.606 fm for Cu nuclei. The 
normalization factor p is set such that J d 3 rpA(f) = A. 
With the above parameters we get po = 0.17 fm -3 . 
The nuclear thickness function 



Ta(x, y) 



dz p A {x,y,z) . 



(10) 



where r = \/ x 2 + y 2 + z 2 , can be used to express both 
the density of wounded nucleons hwn and binary col- 
lisions ubc (see e.g. Ref. Q). Whether the deposited 
energy density or entropy density scales with nwN or ubc 
is not clear from first principles. SPS data suggests that 
the final state particle multiplicity is proportional to the 
number of wounded nucleons. At RHIC energies a viola- 
tion of this scaling was found (the particle production per 
wounded nucleon is a function increasing with centrality. 
This is attributed to a significant contribution from hard 
processes, scaling with the number of binary collisions). 

As in [3[ the shape of the initial energy density distri- 
bution in the transverse plane is parametrized as 

W(x, y, b) = (1 - a) n WN (x, y,b) + a n B c(x, y, b) , (11) 

where a determines the fraction of the contribution from 
binary collisions, and we will study the effect of different 
choices for a. 

For the longitudinal profile we employ the prescription 
used in Refs. pH - ll7j . It is composed of a fiat region 
around ry s = 0, and of half a Gaussian in the forward and 
backward direction: 



exp 



2a 2 



M - TJflat/2) 



The full energy density distribution is then given by 
s(x, y, n.,b) = £ o H(r, s ) W(x, y, b)/W(0, 0, 0) . (13) 



centrality % 


(6) [fm] 


0-5 


1.61 


0-6 


1.76 


0-10 


2.28 


5-10 


2.95 


6-15 


3.47 


10-15 


3.83 


10-20 


4.16 


15-20 


4.55 


15-25 


4.82 


20-30 


5.4 


25-35 


5.92 


30-40 


6.4 



TABLE I. Average impact parameters in selected centrality 
classes from the optical Glauber model. 



The parameters 7/fl a t and a n are tuned to data and will 
be quoted below. 

The impact parameter b is taken to be the average 
impact parameter (b) for a given centrality class and de- 
termined using the optical Glauber model (see e.g. [3). 
The obtained values are listed in Table Q] 

The spectrum of produced hadrons of species i with 
degeneracy gi is given by the Cooper- Frye formula (l9| : 



E 



dN 



dN 



d 3 p dypfdpxc 



•'!> 



g l / f(u^) P »d 3 ^ , (14) 



with the distribution function 



(2tt) 3 exp((u"p M - ih)/Tfo) ± 1 ' 



(15) 



where Tpo is the freeze-out temperature, vf 1 is the flow 
velocity of the fluid, and the energy is given by the in- 
variant expression E = u^p^. In practice, we determine 
the freeze-out hyper-surface geometrically using a trian- 
gulation of the three dimensional hyper-surface that is 
embedded in four-dimensional space-time. The details of 
this algorithm are described in Ref. Q. 



III. HARD PART - MONTE-CARLO JET 
EVOLUTION 

The hard part of the system is described using the 
Monte-Carlo simulation martini Q. At the core of 
martini lies the McGill-AMY formalism for jet evolu- 
tion in a dynamical thermal medium. This evolution 
is governed by a set of coupled Fokker-Planck-type rate 
equations of the form 



(12) dP{p) 



dt 



dk \ P{p+k) 



dT(p+k, k) 
dk 



where dT(p,k)/dk is the transition rate for processes 
where partons of energy p lose energy k. 
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In the finite temperature field theory approach of 
AMY, radiative transition rates can be calculated by 
means of integral equations @ which correctly reproduce 
both the Bethe-Heitler and the LPM results in their re- 
spective limits [20| , Transition rates for the processes 
9 ~^ 99> <?(<?) — ► q{(j)9i an d 9 — > qq are included as well 
as elastic processes employing the transition rates com- 
puted in |7j. Furthermore, gluon-quark and quark-gluon 
conversion due to Compton and annihilation processes, 
as well as the QED processes of photon radiation q — > 
and jet-photon conversions pll |22| are implemented. 

In martini, Eq. is solved using Monte Carlo meth- 
ods, keeping track of each individual parton, rather than 
the probability distributions P. This method provides 
information on the full microscopic event configuration 
in the high momentum regime, including correlations, 
which allows for a very detailed analysis and offers a di- 
rect interface between theory and experiment. 

In a full heavy-ion event, the number of individual 
nucleon-nucleon collisions that produce partons with a 
certain minimal transverse momentum p™ m is deter- 
mined from the total inelastic cross-section, provided by 
pythia, which is also used to generate those individual 
hard collisions. The initial transverse positions of the 
hard collisions are determined by the initial jet density 
distribution 



T A (r ± +b/2)r B (r ± -b/2) 
T AB (b) 



(17) 



which is determined by the nuclear thickness (|10[) and 
overlap functions 



T AB {b) = J d 2 r ± T A (r ± )T B (r ± + b) 



(18) 



The initial parton distribution functions can be selected 
with the hel p o f the Les Houches Accord PDF Interface 
(LHAPDF) [23]. Isospin symmetry is assumed and nu- 
clear effects on the parton distribution functions are in- 
cluded using the EKS98 [24[ parametrization. 

As discussed above, the soft medium is described by 
hydrodynamics. Before the hydrodynamic evolution be- 
gins (r < To), the partons shower as in vacuum. During 
the medium evolution, individual partons move through 
the background according to their velocity. Probabilities 
to undergo an interaction are determined in the local 
fluid cell rest frame using the transition rates and the 
local temperature. The boost into the local fluid rest- 
frame is performed using the flow velocities provided by 
the hydro calculation. 

If a process occurs, the radiated or transferred energy 
is sampled from the transition rate of that process. In 
case of an elastic process, the transferred transverse mo- 
mentum is also sampled, while for radiative processes 
collinear emission is assumed. 

Radiated partons are further evolved if their momen- 
tum is above a certain threshold p m in — 2 — 3 GcV. 
The overall evolution of a parton stops once its energy 
in a fluid cell's rest frame falls below the limit of 4T, 
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PHENIX Cu+Cu 200GeV 0-5%, 5-10%/10. 10-15%/100, 15-207J500, 20-30%/4000 
MUSIC Cu+Cu a=0.1 , b=1 .61 fm, 2.95fm, 3.83fm, 4.55fm, 5.4fm ~~ 
MUSIC Cu+Cu a=1 , b=1 .61 fm, 2.95fm, 3.83fm, 4.55fm, 5.4fm 




FIG. 1. Transverse momentum spectra of negative pions for 
different centrality classes compared to preliminary PHENIX 
data. The collisions go from central (top) to more peripheral 
(bottom). 



where T is the local temperature. For partons that stay 
above that threshold, the evolution ends once they en- 
ter the hadronic phase of the background medium. In 
the mixed phase or a crossover region, processes occur 
only for the QGP fraction. When all partons have left 
the QGP phase, hadronization is performed by PYTHIA, 
to which the complete information on all final partons is 
passed. Because pythia uses the Lund string fragmen- 
tation model [lH, [26| , it is essential to keep track of all 
color strings during the in-medium evolution. For more 
information on the details of martini please refer to Q ■ 



IV. RESULTS 

We begin by presenting results for the soft transverse 
momentum region obtained entirely by MUSIC. We al- 
ways compare three different paramctrizations, which 
essentially differ by the contribution of binary collision 
scaling of the initial energy density distribution (see Eq. 
(|TTT) ) . The details of the parameter sets are listed in Ta- 
ble M 



Fig. [T] shows the transverse momentum spectra of neg- 
ative pions for different centralities compared to prelim- 
inary data from PHENIX (see e.g. [27j]). For the most 
central collisions both paramctrizations describe the ex- 
perimental data reasonably well for pr < 2 GeV. The 
larger the impact parameter, the lower the pr at which 
the description using hydrodynamics fails. For 20-30% 
central collisions the calculation begins to deviate from 
the experimental data at approximately px si 1.5 GeV. 
Generally, the calculation with larger a leads to harder 
spectra. For anti-protons we find a very similar behav- 
ior, shown in Fig. [5] In this case the result is generally 
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a 


to [fm] 


£ [GeV/fm 3 ] 


£Fo[GeV/fm :i ] 


Tpo[MeV] 






0.1 


0.55 


20 


0.14 


w 140 


4 


1.2 


0.5 


0.55 


26 


0.14 


« 140 


4 


1.2 


1 


0.55 


29 


0.14 


w 140 


4 


1.2 



TABLE II. Parameter sets. 



PHENIX Cu+Cu 200GeVpO-5%,5-10%/10, 10-15%/100, 15-20%/500, 20-30%/4000 *&* 
MUSIC Cu+Cu R=0.1 , b=1 .76fm, 2.95fm, 3.83fm, 4.55tm, 5.4fm ~~ 
MUSIC Cu+Cu a=1 , b=1 .76 fm, 2.95frru 3.83fm, 4.55fm, 5.4fm 




300 
250 
200 



5^ 



PHOBOS prelim. 0-6%, 6-15%, 15-25%, 25-35% 
MUSIC Cu+Cu o=0.1, b=1.76fm, 3.47fm, 4.82fm, 5.92fm 
MUSIC Cu+Cu a=0.5, b=1.76fm, 3.47fm, 4.82fm, 5.92fm 
MUSIC Cu+Cu a=1 , b=1 .76fm, 3.47fm, 4.82fm, 5.92fm 




FIG. 2. Transverse momentum spectra of anti-protons for 
different centrality classes compared to preliminary PHENIX 
data. The experimental data is underestimated because of 
the assumption of chemical equilibrium. 



FIG. 3. Centrality dependence of pseudorapidity distributions 
of charged hadrons compared to preliminary PHOBOS data 
[3uj ] (open symbols). Calculations using a = 0.1 describe all 
centrality classes very well, while larger a lead to deviations 
for larger centralities. 



below the data because of the assumed chemical equilib- 
rium (see Ref. Q), the overall shape of the spectra is 
however well reproduced for low px- 

Although we find harder spectra in the case of a = 1 
it is unclear which parametrization is favored solely by 
studying the pr spectra of pions and anti-protons. The 
centrality dependence of charged particle pseudo-rapidity 
distributions in Fig. [3] reveals a much more obvious dif- 
ference. While in the case a = 0.1 the experimental data 
is very well described for all shown centrality classes, the 
calculations using a = 0.5 and a = 1 only describe the 
most central data to which the parameters were adjusted, 
but differ increasingly with increasing impact parameter. 

It is remarkable how well the px integrated pseudo- 
rapidity distributions are described even for very large 
pscudorapidities and impact parameters. This indicates 
that for low pt, i.e., pr ^ 1 GeV, hydrodynamics works 
well, hinting at that the low momentum modes equili- 
brate fast enough even in small systems, while higher 
momentum modes do not. 

Next, we present the elliptic flow parameter vi and 
compare to STAR data j28[. Because we do not include 
fluctuations, we do not expect the centrality dependence 
of elliptic flow to be reproduced. See Ref. [29| for a 
demonstration of the importance of fluctuations on el- 
liptic flow observables. In Fig. @] we show results for 
10-20% central events, a centrality class in which V2(pt) 
was found to be well described by smooth initial condi- 




FIG. 4. Elliptic flow coefficient vi as a function of transverse 
momentum for 10-20% central events compared to STAR data 

M. 



tions, at least for Au+Au collisions [29(. The dependence 
on the a parameter is as expected. Larger a meaning 
higher initial eccentricity, leads to stronger elliptic flow. 
To draw further conclusions from this comparison a de- 
tailed event- by- event study is needed. 

Extracting the evolution information from the hydro- 
dynamic simulation, i.e., the temperature and flow veloc- 
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FIG. 5. Transverse momentum spectrum of negative pions 
from the hydrodynamic calculation MUSIC (dashed line) and 
martini using MUSIC input (solid line) compared to data from 
PHENIX (preliminary) and STAR HJ. 



ities as functions of time and spatial coordinates, we can 
use this information as an input for the simulation of the 
evolution of the hard degrees of freedom, martini. 

This way we can compute the hard part of the spec- 
trum and combine results for the soft and hard part as 
done in Fig.[S] for negative pions. The additional free 
parameter that regulates the amount of energy loss in 
martini is the strong coupling constant, which we fix 
at a s = 0.32. The value of transverse momentum up 
to where hydrodynamics alone provides a good descrip- 
tion and beyond which the contribution from jets be- 
comes important decreases with decreasing system size. 
For central Au+Au collisions the value lies between 3 
and 4 GeV Q, while Fig. [5] shows that for Cu+Cu it lies 
around 2 GeV for central collisions, and from Fig.Q] one 
can read off that it decreases more when going to larger 
centralities. 

Next, we study the nuclear modification factor for neu- 
tral pions, defined by 



Raa 



1 dN AA (b)/d 2 PT dy 
iVcou(6) dN pp /d 2 p T dy 



(19) 



where N co u(b) is the number of binary collisions at 
given impact parameter b. The reference spectrum 
dN pp /d 2 pTdy for p+p collisions was computed in Q. 

Figures [5] to [5] show the nuclear modification factor 
Raa for different centrality classes. The value of a s was 
chosen to match the experimental data in the most cen- 
tral collisions best, but when going to larger centrality 
classes the impact parameter is the only quantity that 
is changed. We see that all hydrodynamic backgrounds 
(a e {0.1,0.5,1}) lead to good agreement with the ex- 
perimental data in the two most central bins (Figs. [6] and 
[7]). All calculations, but most significantly the one using 
a = 0.1, deviate from the 20-30% and 30-40% central 
experimental data. This indicates that the combination 
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FIG. 6. Nuclear modification factor for neutral pions in 0- 
10% central Cu+Cu collisions compared to preliminary data 
from PHENIX. 
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PHENIX prelim. Cu+Cu 200 GeV, 10-20% central 
AMY + el 3+1 D MUSIC Hydro a s =0.32a=0.1 b=4.16fm 
AMY + el 3+1 D MUSIC Hydro a s =0.32 a=0.5 b=4.1 6 fm 
AMY + el 3+1 D MUSIC Hydro <x s =0.32 a»1 .0 b=4.1 6 fm 
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FIG. 7. Nuclear modification factor for neutral pions in 10- 
20% central Cu+Cu collisions compared to preliminary data 
from PHENIX. 



of ideal hydrodynamics with our energy loss mechanism 
begins to fail to describe Cu+Cu systems at centralities 
> 20%. Remember that hydrodynamics using a = 0.1 
described the bulk data well, while deviations from ex- 
perimental data at more peripheral collisions are large 
for a = 0.5 and a = 1 (Fig.©. 

On the other hand a = 1 works best for Raa- This 
indicates that the energy loss mechanism itself leads to a 
too weak centrality dependence when using the "correct" 
hydrodynamic background (a = 0.1). Using a hydro- 
background with a larger centrality dependence of the 
density (a = 1) then works against this problem. Note, 
however, that for all a the results for Raa lie within the 
experimental error bars. The effect of non-zero viscosity 
and event-by-event fluctuations [20l ] remains to be stud- 
ied. 

Finally we study high momentum photon production 
in Cu+Cu collisions. Photons in the high px region pro- 
duced in nuclear collisions are dominantly direct photons, 
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PHENIX prelim. Cu+Cu 200 GeV, 20-30% central 
AMY + el 3+1 D MUSIC Hydro <x s =0.32 <x=0.1 b=5.4 fm 
AMY + el 3+1 D MUSIC Hydro <x s =0.32 a=0.5 b=5.4 fm 
AMY + el 3+1 D MUSIC Hydro a s =0.32 a=1 .0 b=5.4 fm 
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PHENIX prelim. Cu+Cu 200 GeV, 0-10% central 
AMY + elastic + MUSIC Hydro, a s =0.32, b=2.28 fm, a=1 
AMY + elastic + MUSIC Hydro, a s =0.32, b=2.28 fm, a=0.1 
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FIG. 8. Nuclear modification factor for neutral pions in 20- 
30% central Cu+Cu collisions compared to preliminary data 
from PHENIX. 



FIG. 10. Nuclear modification factor for photons in 0-10% 
central Cu+Cu collisions compared to preliminary data from 
PHENIX. 
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PHENIX prelim. Cu+Cu 200 GeV, 30-40% central 
AMY + el 3+1 D MUSIC Hydro ct s =0.32 ct=0.1 b=6.4 fm 
AMY + el 3+1 D MUSIC Hydro ct s =0.32 ct=0.5 b=6.4 fm 
AMY + el 3+1 D MUSIC Hydro a s =0.32 a=1 .0 b=6.4 fm 
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FIG. 9. Nuclear modification factor for neutral pions in 30- 
40% central Cu+Cu collisions compared to preliminary data 
from PHENIX. 



fragmentation photons, and jet-plasma photons. Direct 
photons are included in PYTHIA. Apart from leading or- 
der direct photons, pythia produces additional photons 
emitted during the vacuum showers, which in part over- 
lap with photons from next-to-leading order calculations. 
Fragmentation photons are also part of those produced 
in the shower. For heavy-ion reactions martini adds the 
very relevant jet-mcdium photons from photon radiation 
q — > ^7 (bremsstrahlung photons) and jet-photon con- 
version via Compton and annihilation processes. The full 
vacuum shower is included in the calculation - most of the 
shower photons will be emitted before the medium has 
formed and the parton has realized that it has formed. 
For more details on the implementation see [||. 

In Fig. [10] we present the nuclear modification factor 
for photons compared to preliminary data from PHENIX 
(also see [27|) in central collisions. Again, the reference 
spectrum from p+p collisions was computed in Q. The 



data is consistent with R\ A = 1 and so is the calculation 
for a wide range of px- The rise above one for px ^ 
6 GeV stems from the contribution of in-medium photon 
production. Unfortunately the current experimental data 
does not allow to distinguish whether this contribution 
is present or not. 



V. CONCLUSIONS AND OUTLOOK 

We studied the applicability of both a 3+1D hydrody- 
namic model and the McGill-AMY jet energy loss and 
evolution scheme to Cu+Cu collisions of different ccn- 
tralitics. Because the created system is smaller than in 
Au+Au collisions the limits of the models are already 
reached at a smaller impact parameter. The precise de- 
termination of these limits and the determination of ob- 
servables that are sensitive to the failing of the models 
was the main objective of this work. 

The hydrodynamic evolution was computed using 
the relativistic 3+1D implementation of the Kurganov- 
Tadmor scheme, dubbed MUSIC, while the evolution of 
the hard degrees of freedom was simulated using the 
Monte-Carlo simulation martini. With these, we were 
able to compare calculations to experimental observables 
over the full range of measured transverse momentum. 

As expected, the quality of the description of trans- 
verse momentum spectra of identified hadrons using only 
hydrodynamics depends on the system size. While in 
Au+Au collisions the spectra in central collisions are typ- 
ically well described up to px ~ 3 GeV, in Cu+Cu we find 
deviations at approximately 2 GeV for central collisions, 
and lower px for larger centralities. 

The hydrodynamic calculation of the centrality depen- 
dence of pscudorapidity distributions agrees very well 
with the experimental data when we use 90% wounded 
nucleon scaling and 10% binary collisions scaling of the 
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initial energy density distribution. The parameter set 
using 100% binary collision scaling does not reproduce 
the centrality classes above 10%. We therefore conclude 
that mostly wounded nucleon scaling provides a better 
description of experimental data. In this case, the result 
shows that the low transverse momentum modes that 
dominate the distribution arc thermalized, even in 30- 
40% central collisions, which is in agreement with the 
result for the transverse momentum spectra. 

Including the contribution from jets by employing 
martini, which uses the hydrodynamic evolution as in- 
put, we were able to describe the complete pion spec- 
tra in central collisions. The study of Raa for pions 
allows to pin down the point where the combined ideal 
hydrodynamics+McGill-AMY evolution breaks down. 

The trend seen in the description of the centrality de- 
pendence of Raa agrees least with experimental data 
when using a = 0.1, which is preferred by the low mo- 
mentum spectra. While calculations for 0-10% and 10- 



20% central bins agree well with experimental data, devi- 
ations begin for centralities larger than 20%. This might 
indicate that the used energy loss mechanism does not 
have the correct density and/or system size dependence 
for the small systems studied here. The current error bars 
do not permit a more precise quantitative statement. 
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